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. Abstract. Given the increasing availability of lattice data for (unquenched) QCD with Nf = 2, it 

£jfy is worth while to check whether the generated vacuum significantly deviates from the quenched 

one. I discuss a specific attempt to do this on the basis of topological susceptibility data gained at 
various sea-quark masses, since for this observable detailed predictions are available. The upshot is 
that either discretization effects in dynamical simulations are still untolerably large or the vacuum 
C*~) \ structure in 2-fiavour QCD substantially deviates from that in the theory with 3 (or 2+1) light quarks. 



> : INTRODUCTION 

m 

O ■ In a historical perspective, the path towards phenomenological predictions of QCD by 

means of lattice techniques involves three steps: Pure Yang-Mills theory (where there is 
just glueball physics), quenched QCD (where the vacuum is the one in the SU(3) theory, 
but observables may involve so-called current-quarks) and full QCD (where the fermion 
determinant with the dynamical sea-quarks accounts for the quark loops in the vacuum). 
Today, the lattice community makes the final push towards full QCD, despite the fact 
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that state-of-the-art simulations are modestly announced as "partially quenched" which 
means that the sea- and current-quark masses in the (euclidean) generating functional 

Z[fi,r|] = [DA e~ s ° n det (#+ m sea) e*p( / fj— - — r|) (1) 



are (in general) unequal and in most cases significantly heavier than the physical u- and 
J-quarks, so that phenomenological statements require a twofold extrapolation. 

Since the finite sea-quark mass constitutes the key ingredient in this ultimate step 
(note that the determinant turns into a constant for m sea — > °°, hence (1) reduces to 
the quenched generating functional in that limit), an obvious task is to check whether 
these "partially quenched" or "full" QCD simulations exhibit the change in the vacuum 
structure expected to occur if the fermions are "active" (i.e. if the back-reaction of 
the "dynamical" fermions on the gauge background is taken into account). The prime 
observable used to distinguish the respective vacua is the topological susceptibility 

X(m sea ) = , (2) 

with q the (global) topological charge, because detailed theoretical predictions show 
that % behaves rather different in the quenched (m sea — > °°) and chiral (m sea — > 0) limits, 



respectively. Even though in the lattice-regulated theory (and with certain definitions of 
the topological charge operator) q may be somewhat ambiguous on the level of a single 
configuration, the moment of the ^-distribution which enters (2) can be measured with 
controlled error-bars, and as a purely gluonic object the resulting % = %(m sea ) encodes 
nothing but the vacuum structure of the theory. 

Below, I give a quick survey of recent lattice determinations of % at various sea-quark 
masses in N/ = 2 QCD, I discuss the available continuum knowledge of the functional 
form % = %(w S ea), and I present a non-standard lattice determination of the quenched 
topological susceptibility and the chiral condensate E based on it. The outcome is 
that either certain observables in todays phenomenological studies with light dynamical 
quarks suffer from large discretization effects or - the more speculative view - that the 
low-energy structure of QCD with Nf = 2 is substantially different from that with N/ = 3. 

LATTICE DATA 

I start with a quick survey of recent lattice data for the topological susceptibility in QCD 
with 2 dynamical flavours; the selection reflects nothing but my personal awareness. 

CP-PACS: The CP-PACS collaboration has simulated full QCD on several grids 
at various (|3, K)-values, using an RG-improved gauge action and an 0(a) -improved 
fermion action with mean-field values for the associate csw coefficients. Below, I con- 
centrate on the data generated on the 24 3 x 48 lattice at [3 = 2.1 with LW-cooling [1]. 

UKQCD: The UKQCD collaboration has simulated full QCD on a 16 3 x 32 grid at 
various (|3,K)-values, using the standard (Wilson) gauge action and an 0(a) -improved 
fermion action with the non-perturbative values for the associate csw coefficients [2]. 

SESAM/TXL: The SESAM/TXL collaboration has simulated full QCD on two grids 
(16 3 x 32 and 24 3 x 40) at several ((3, K)-values, combining the unimproved (Wilson) 
gauge action with unimproved (Wilson) fermions (i.e. setting csw = 0) [3]. 

Thin link staggered: Trusting a continuum identity for the relationship between the 
2- and the 4-flavour functional determinant, the staggered fermion action may be used 
to simulate QCD with N/ = 2. There are data by the Pisa group [4] and by A. Hasenfratz 
based on configurations by the MILC collaboration and the Columbia/BNL project [5]. 

Fig. 1 displays the data, along with continuum constraints to be discussed next. 



CONTINUUM KNOWLEDGE 

As mentioned in the introduction, the data for the topological susceptibility % versus the 
sea-quark mass m = m sea prove useful to test the vacuum structure, because continuum 
QCD provides us with detailed predictions: There are analytic upper bounds for %(m) 
at both asymptotically small and large sea-quark masses and there is a "semi-analytic" 
formula for %(m) valid at intermediary quark masses (where the bulk of the lattice data 
reside). The only caveat is that these bounds hold true in the continuum limit, but so far 
no continuum extrapolation for %(m) in N/ = 2 QCD is available yet. Before stating the 
complications due to this, the continuum functional forms shall be discussed. 




FIGURE 1. Topological susceptibility versus quark mass (each in dimensionless units) in Nf = 2 QCD 
with Wilson-type (left) or staggered (right) sea-quarks [1, 2, 3, 4, 5]. For comparison, la bands indicating 
the constraints in the deep chiral regime (based on F % =93 ± 1 MeV) and in the heavy sea-quark (quenched) 
limit (from = (200 ± 5 MeV) 4 ) are shown as well as the associate continuum band (5) (full line). 



Asymptotically small sea- quark masses: For m<^ Aqcd and to leading order in the 
chiral expansion the axial WT-identity yields (see Refs. cited in [6] for details) 

X(m) = ^(l+0(m)) = *@£(l+0(m)) = Xo (l+0(m)) , (3) 

where, in the second equality, the Gell-Mann-Oakes-Renner relation has been assumed. 

Asymptotically large sea-quark masses: For m 3> Aqcd the topological suscepti- 
bility gradually approaches its quenched counterpart (see Refs. cited in [6] for details) 

X(m) = 3Coo(l+0(l/m)) = (200 ±5 MeV) 4 , (4) 

Intermediate quark masses: For other quark masses (i.e. for (M 7r ro) 2 G [1 .5, 15] or so, 
with the Sommer scale ro = 0.5 fm throughout) neither one of the asymptotic predictions 
is applicable. Fortunately, there is the "reduced" interpolation formula 

3C(m) = l/(l/xo + l/X~) (5) 

with %o defined in (3), which is, of course, not exact but represents an "educated guess"; 
it follows either from the chiral Lagrangian together with pure entropy considerations 
(which makes it very robust) [6] or from large-iV c - arguments [7]. 
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FIGURE 2. Topological susceptibility versus quark mass in Nf = 2 QCD with Wilson-type (left) or 
staggered (right) sea-quarks [1, 2, 3, 4, 5] together with naive fits of the susceptibility curve (5) to the data, 
neglecting possible discretization effects. The associate values for £2 an d 5C~ are tabulated in Table 1. 



NAIVE EVALUATION OF £ 2 AND x~ 

Disregarding possible lattice artefacts, one may fit the available data to the continuum 
curve (5) and extract E and from the fit parameters [6]. It is worth emphasizing that 
this evaluation of E is distinct from the usual fermionic determination which is via the 
trace of the Green's function of the Dirac operator at various current-quark masses and 
extrapolating (after proper renormalization) to the physical (or chiral) point. Regardless 
how convincing this sounds, the results as tabulated in Table 1 look rather devastating: 
While our value for the quenched topological susceptibility %oo — (200± lOMeV) 4 nicely 
agrees with previous direct determinations in the SU(3) -theory, the suggested value for 
the (full) chiral condensate in the chiral limit E~ (450± 100 MeV) 3 dramatically exceeds 
(by more than a factor 2) the "phenomenological" value E~ (288 MeV) 3 (which follows 
from the GOR-relation with m v (MS,^ = 2GeV) -3.5 MeV [8]). 

Looking back at Fig. 1, one may argue that this hardly comes as a surprise, since 
both the "Wilson-type" and the "staggered" data sets are much more likely to violate 
the linear upper bound in the deep chiral regime than the flat ceiling in the heavy-(sea)- 
quark limit. Besides, Fig. 1 tells us how important it is to compare the data to the right 
prediction: Knowing nothing but the chiral constraint (3), one might be tempted to say 
that the data are in nice agreement with the leading order chiral prediction. The outcome 



TABLE 1. Naive determination of E and %oo from full QCD vacuum data 
with (5), using ro = 0.5 fm and the GOR-relation to convert to physical units. 





CP-PACS 


UKQCD 


SESAM * 


PISA 


BOULDER 


{Fnr ) 2 /4 


0.0417 


0.0216 


0.1600 


0.1728 


0.0441 


F n [MeV] t 


161. 


116. 


316. 


328. 


166. 


zV 3 [MeV] 


415. 


334. 


650. 


667. 


423. 


4 


0.0616 


0.0781 


0.0587 


0.0343 


0.0748 


xi /4 [MeV] 


197. 


209. 


194. 


170. 


206. 



* All data get equal weights, otherwise the best direct fit is almost flat (cf. Fig. 2). 

^ Using the convention in which fj[~93MeV in nature; note that - except for the 
one in the UKQCD column - all entries are substantially larger than that value. 



of our analysis shows that there is absolutely no point in comparing lattice data gained at 
(M n ro) 2 > 1.5 to the leading order chiral behaviour (3), because the "true" prediction (5) 
is substantially lower: If lattice data at (M JI ro) 2 = 2.5 are found to be in "nice agreement" 
with the chiral prediction (3) based on the phenomenological value E~ (288MeV) 3 , then 
it means that they are ~50% in excess of what they should be. 

The bottom line is that todays full QCD simulations (with both Wilson-type and 
staggered sea-quarks) - if taken at face value - do show unquenching effects in their 
vacuum structure but, in general, far less than expected at their respective sea-quark 
masses. Unpleasant as it is, we are invited to think about the reasons for this finding. 

TWO ALTERNATIVES 

There are two main reasons why E as determined via fitting full QCD topological 
susceptibility data to (5) could substantially exceed the standard phenomenological 
value E~ (288MeV) 3 while the simultaneously determined %^ takes a regular value. 

Large lattice artefacts: The simple reason is that discretization effects could be large, 
since for the case of the topological susceptibility finite-volume effects are analytically 
shown to be well under control in most of todays simulations [6]. With the ascent of 
fermion actions which satisfy the Ginsparg-Wilson relation it makes sense to disentangle 
chirality violation effects from ordinary scaling violation effects, and in a recent paper 
A. Hasenfratz has shown some evidence [5] that one should primarily suspect the former 
type of discretization effects to give rise to the excessive £ values listed in Table 1 . 

Proximity of the "conformal window": The more "exotic" view is that the fitted 
value for %o( m ) — Em/A^- is appropriate for Nf=2 QCD, while the standard phenomeno- 
logical evaluation - even if it involves (non- strange) pions only - concerns QCD with 
3 (2+1) light flavours. In this scenario the way % depends on m (m<^AQco) has an 
A^-dependence beyond the one indicated in (3), i.e. the low-energy constant E depends 
(strongly) on N f , e.g. E 2 ~ (450MeV) 3 while E 3 ~ (288 MeV) 3 . The latter gap could then 
be interpreted as a hint that the "conformal window" (where, for appropriate Nf, one has 
E^ <C Aq CD and chiral symmetry is primarily broken through higher dimensional con- 
densates) might be "close" - see [6] for a discussion and some references. 




FIGURE 3. Schematic representation how discretization effects typically affect topological suscepti- 
bility data in full QCD with Wilson-type (left) and staggered (right) sea-quarks, if simulations are run at 
fixed P and fixed m, respectively. In either case discretization effects tend to enhance the effective F % or £ 
(in particular for large lattice spacing a) and they reduce the associate %oo by an amount oca 2 . The lattice 
susceptibility (6) is the reduced mean of the modified functions (fat dots) and may, for 2 < (M^ro) 2 < 8, 
show little variation with m and lie substantially above the expected continuum curve (black line). 

REFINED EVALUATION OF £ 2 AND %«> 

In the following, I concentrate on the first alternative and show that - in the absence of 
data suitable for a continuum extrapolation - basic knowledge regarding the dominant 
discretization effects allows for a more sophisticated evaluation of the parameters in (5). 

The key observation on which this analysis relies is that the leading lattice artefacts 
in both elements of (5) - the chiral piece (3) and the quenched piece (4) - are known: 
On the chiral side the dominant effect is chirality violation (for a discussion see [5, 2]), 
i.e. F n ro = F n ro(\ + const (a/ro) p ) with p reflecting the fermion formulation. On the 
quenched side scaling violations are known to result in = X°° r o — 0.208(a/ro) 2 
(with a known coefficient !) [2, 1]. Combining all the ingredients, one ends up with 

P 4 = \/{2N f /[(M n ro) 2 (F n ro) 2 (\ + const(a/ro) p ) q } + l/[X^-0.208(a/r ) 2 ] } , (6) 

where q may be chosen between 1 and 2, since (l + 0(a p )) 2 = \-\-0(a p ); I use q = 2. 
A qualitative picture how these modifications affect the measured topological suscepti- 
bility is drawn in Fig. 3. In this respect it is important to know that in a series of full QCD 




FIGURE 4. Topological susceptibility versus quark mass in Nf = 2 QCD with Wilson-type (left) or 
staggered (right) sea-quarks [1, 2, 3, 4, 5] together with fits of the functional form (6) to the data (dotted 
lines). This time, the values for £2 an d 5C°° suggest "reasonable" continuum curves (full lines, c.f. Table 2). 

simulations at fixed [3 the lattice spacing will shrink if the (sea-)quark mass gets reduced 
(which is often the case in studies with Wilson-type sea-quarks; the only exception is 
the one by UKQCD, where [3 gets relaxed when K is increased in such a way that a stays 
approximately constant), whereas if one works in a staggered formulation at fixed quark 
mass m the lattice will get coarser as one approaches the chiral limit - eventually the 
measure resembles more that of a theory with a single pseudo-Goldstone rather than one 
with Nj — 1 pion type degrees of freedom as in the continuum [5]. 

In order to make use of this knowledge (i.e. to utilize (5) to fit the data) one has to 
decide on the parameters (p, const) showing up in (6). The former choice is relatively 
easy - 1 take p cp-pacs = 1 • 5 , p tjkqcd = 2 , p ses am = 1 , P pis a = P boulder = 2 to account 
for the formulation and the different strategies regarding csw- The latter choice - which 
value "const" shall be given - is more delicate: Ideally, one would like to determine it 
directly from the data. However, it turns out that the quality of the data at hand is not suf- 
ficient to allow for an additional (i.e. third) fitting parameter. A reasonable option would 
be to determine it from conventional F n measurements on the individual ensembles. A 
simpler option is to use (6) twice - in a first round "const" is given a likely value by fit- 
ting it while ((F n r ) 2 /4, X°°4) is held fixed at (0-014, 0.066); in a second round the latter 
get adjusted while "const" is kept fixed at the previously determined value. Obviously, 
with this simpler option the final outcome for ((F % ro) 2 / '4,%^^) reflects, to some extent, 



TABLE 2. Refined determination of E and from full QCD vacuum data 
with (6), using ro =0.5 fm and the GOR-relation to convert to physical units. 
For a cautionary statement regarding the fitting procedure see text. 





CP-PACS 


UKQCD 


SESAM 


PISA 


BOULDER 


(^r ) 2 /4 


0.0174 


0.0134 


0.0299 


0.0340 


0.0106 




104. 


91. 


137. 


146. 


81. 


E^pVleV] 


311. 


284. 


372. 


388. 


263. 




0.0737 


0.0845 


0.0661 


0.0562 


0.1466 


xi /4 [MeV] 


206. 


213. 


200. 


192. 


244. 



the corresponding initial values. The simplest option is just to set "const" to a generic 
value like 1. My person choice is to take the arithmetic average of the "const" values 
suggested by the last two options and to use that value to fit for (F n ro) 2 /A and X^r^. 

The result of this exercise is shown in Table 2 and Fig. 4, where dotted lines represent 
the lattice curve (6) with (F n ro) 2 /4 and %^ r o adjusted such as to make it go through 
the data points while full lines indicate the associate continuum curve (5). The data in 
Table 2 should be taken with a grain of salt since, as explained above, there is a remnant 
trace of the initial values in the final fitting parameters (F n ro) 2 /4 and %^ r o- Nontheless, 
the result is interesting because it supports the standard view that QCD with Nf = 2 
is not in the "conformal window" and that its low-energy structure agrees with that 
suggested by phenomenological investigations in "real" (2+1 flavour) QCD [9], i.e. at 
least for Nf = 2 chiral symmetry is predominantly broken through a distinctively nonzero 
condensate (for references to an alternative scenario see [6]) and E2 as suggested by 
Table 2 is compatible with the value from the GOR-relation, E2+1 ~ (288MeV) 3 . 

From a lattice perspective it reassuring to see that through a simple ansatz for the 
dominant discretization effects todays state-of-the-art simulations (which, from a naive 
perspective, seemed to support an almost flat topological susceptibility curve and hence 
to reproduce - in spite of all unquenching efforts - a more or less quenched vacuum 
structure) may, in fact, be shown to give supportive evidence in favour of a decreased 
topological susceptibility near the chiral limit and hence "bear" the knowledge of the 
difference between the quenched and the unquenched vacuum structure in them. The 
ultimate goal is, of course, to make discretization effects sufficiently small so that the 
expected continuum pattern of the vacuum structure gets visible in the raw data already. 
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